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Two-Loop Effects and Current Status of 
the 4He+ Lamb Shift 



U. D. Jentschura and M. Haas 



Abstract: We report on recent progress in the treatment of two-loop binding corrections to 
the Lamb shift, with a special emphasis on S and P states. We use these and other results in 
order to infer an updated theoretical value of the Lamb shift in ^He^. 

PACS Nos.: 12.20.Ds, 31.30.Jv, 06.20.Jr, 31.15.-p 

Resume : Nous examinons le progres recent concemant le traitement des corrections a I'ordre 
eleve des diagrammes a deux boucles contribuants aux deplacement de Lamb, specialement 
en ce qui concerne les etats S et P. Par consequent, on deduit de nouvelles precises valeurs 
theoriques pour le deplacement de Lamb en ^He^. 

1. Introduction 

Recently, the higher-order two-loop corrections to the Lamb shift have been studied rather inten- 
sively, both within the Za-expansion (see [1-8] and references therein) as well as within the nonpertur- 
bative (in Za) numerical approach, as described in Refs. [9-15]. In the current note, we review some 
recent progress for the so-called Bqq coefficient, which is generated by the entire gauge-invariant set 
of two-loop diagrams, as depicted in Fig. 1. We also review some very recent progress [16] regarding 
the i?6o coefficient for general excited hydrogenic states with nonvanishing angular momentum, with 
a special emphasis on states with P symmetry. 

Applications of the recent progress to high-precision spectroscopy are numerous. As one example 
of current interest, the status of the "^He^ Lamb shift (IS and 25" states) is summarized, based on 
the recent analytic and numerical results, and on information about further known contributions to the 
Lamb shift from the literature (see in particular [7, 17, 18] and references therein). 

2. Two-Loop Results 

The two-loop energy shift of an atomic level in a hydrogenlike atomic system reads (in units with 

^ = c = eo = 1) 



where ttIc is the electron mass, and H isa dimensionless function. In the current Section of this article, 
we are primarily concerned with recently obtained [8, 16] results for the normalized (or "weighted") 
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Fig. 1. Fcynman diagrams for the two-loop self-energy corrections, separated into subsets i-iv according to 
Ref. [8]. Subset i is the pure two-loop self-energy, subset ii comprises the vacuum-polarization insertion into 
the virtual-photon line of the one-photon self-energy, subset Hi contains vacuum-polarization corrections to the 
electron line in the one-photon self-energy, and subset iv contains remaining vacuum-polarization effects. 



difference H{nS, Za) — -ff (15, Za) of S states, whose importance for the determination of funda- 
mental constants has been stressed in Refs. [19-21], and for individual P states. 

For these states and/or combinations of states, the first nonvanishing terms in the semi-analytic 
expansion of the dimensionless function H{Za) in powers for Za and In(Za) read as follows. 



H{Za) = B40 + {Zaf\ Be2 ln^[{Za)-^] + Bei ln[(Za)-2] + B, 



60 



(2) 



The first index of the B coefficients marks the power of Za, whereas the second corresponds to the 
power of the logarithm ln[(ZQ!)~^]. For individual S states, we only mention here the existence of a 
B50 coefficient [1-3,5], which goes beyond the coefficients listed in (2). 
For the normahzed difference of S states, we have [6] 



Be2inS) - B62(1S) = ^f^ + 2_.l+^- ln(n) + *(n; 



(3) 



where 7 = 0.577216 ... is Euler's constant, and ^ is the logarithmic derivative of the Ganrnia function. 
The normalized difference for Bqi reads [6] 



Bei {nS) - Bei (15) = - [N{nS) - N{1S)] 



+ 



304 32 , 

ln(2) 

135 9 ^ ^ 



1 



1 



n 

The normalized difference of the nonlogarithmic term can be expressed as [8] 

BeoinS) - Beo{lS) = bUnS) - 6^(15) + A(n), 



7 - ln(n)-|-*(n) 



(4) 



(5) 
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where A{n) is an additional contribution beyond the n-dependence of the two-loop Bethe logarithm 
6i. The result for A is [8, 22], 

./ X /^38 4, r,.. ox r.M 337043 94261 902609 

A{n) = — - - ln(2) [N{nS) - N{1S)] 



V45 3 ' 7 ^ ^ ^ ^ ^' 129600 21600n 129600n2 
I - 16 4 \ , 2,„x / 76 304 76 \ , / 53 35 419 \ 

+ ( o - 9^ + g^^j 1- (2) + (-45 + 135^ - 135^j 1^(2) + (-15 



3 



X 



+ (iSS + ¥ ^^'^^^ " S ^''^^^ + ^^^^^ " y ^^^^ " ^^^^0 ^ ^^""^ ~ ^''^''^^ ■ 

Here, A''(n) is a nonlogarithmic term generated by a Dirac-(5 correction to a one-loop Bethe logarithm, 
as calculated in Ref. [23]. Of course, ((s) = X^J^i '^ ^ the Riemann zeta function. 
For P states, we have the known results [4, 8, 22] 

Be2{nP) ^ " ^ 



27 n2 



B.i(nP„.) = 3 N(nP) + _ - - l„2j . (7) 

The results for the nonlogarithmic terms of P can be inferred on the basis of Eq. (8.1) of Ref. [8] and 
the two-loop Bethe logarithms for P states (see [16] and Table 2), 



B6o{nPi/2) = bL{nP) + /34(riPi/2) + ^^{nPi/^) 



38 4, 

ln(2 

45 3 ^ 



27517 209 

N(nP) 

^ ' 25920 288n 



1223 4 n2-l,2,„, 38n2-l, /25 3 9 , , 

/ 9151 1 1009 \ , / 25 3 9 \ , 
+ -10800 - Tn + 1200^? ^^'^ + "24 " 8^^ + 8^? ^^^^ ' ^^^^ 



B6o{nP3/2) = bL{nP) + p4{nP3/2) + p5{nP3/2) + 



38 4, 

45-3^^^^) 



73321 185 
^ 103680 1152n 



8111 4n2-l 2/„x lln^-l /299 3 53 

+ 25920^? + 27^^^ - 81^ ^^(2) + " 8^^ " 20^? > ^^^^^^^^^ 

/ 24377 1 3187 \ / 299 3 53 \ 

+ (-21600 + T6;i - 3600;?j ^^'^ + (" 320 + ^ + 80^? j ^^'^ " ^'""^ 

In these formulas, /34 and 05 are low-energy spin-dependent contributions, defined in Eq. (4.21) of 
Ref. [8], whose numerical values may be inferred from one-loop calculations [23, 24]. 

The evaluation of the two-loop Bethe logarithmfor 15* and 25* has been discussed in Ref. [25], and 
for 35-65 in Ref. [26]. For IS" and 2S, there is no ambiguity in the definition of the Bethe logarithm. 
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Table 1. Total numerical values of the two-loop Bethe logarithms 6i, for 5 and P states, broken 
down for the principal- value contribution 6i and the squared-decay term d^Beo- 



level 


bL 




bL 


level 


bL 


S'^Beo 


bL 


15 


-81.4(3) 


0.0 


-81.4(3) 










25 


-66.6(3) 


0.0 


-66.6(3) 


2P 


-2.2(3) 


-0.008 


-2.2(3) 


35 


-63.5(6) 


-0.071 


-63.6(6) 


3P 


-2.5(3) 


-0.177 


-2.7(3) 


45 


-61.8(8) 


-0.109 


-61.9(8) 


4P 


-2.8(3) 


-0.243 


-3.0(3) 


55 


-60.6(8) 


-0.129 


-60.7(8) 


5P 


-2.8(3) 


-0.276 


-3.1(3) 


65 


-59.8(8) 


-0.141 


-59.9(8) 


6P 


-2.9(3) 


-0.295 


-3.2(3) 



Table 2. Numerical values for the weighted difference of Bqo for 5 states, and for individual 
P states. 



level Seo (n5) -S6o( 15) level B6o(2Pi/2) level B6o(2Pv2) 



25 


15.1(4) 


2Pl/2 


-1.6(3) 


2P3/2 


-1.8(3) 


35 


18.3(7) 


3Pl/2 


-2.0(3) 


3P3/2 


-2.2(3) 


45 


20.0(10) 


4Pl/2 


-2.4(3) 


4P3/2 


-2.5(3) 


55 


21.2(11) 


5Pl/2 


-2.4(3) 


5Pv2 


-2.5(3) 


65 


22.0(11) 


6Pl/2 


-2.5(3) 


6P3/2 


-2.6(3) 



which can roughly be explained as follows: essentially, the two-loop Bethe logarithm results from a 
renormahzed integration over two photon energies. Both of these integrations are free of singularities 
for 15 and 25. However, for all higher excited 5 states and all P states, one incurs real (rather than 
imaginary) contributions to the energy shift from the product of imaginary contributions due to singu- 
larities along both photon integrations (these are "squared decay rates" in the sense of Ref. [27]). It is 
thus necessary to make a clear distinction between the singularity-free, principal- value part bL and a 
real part S'^B^o, which is incurred by "squared" (or, more precisely, products of) imaginary contribu- 
tions from the pole terms. We write 

bL = bL + S^Bqo , (9) 

where 6l is obtained as the nonlogarithmic energy shift stemming from the nonrelativistic self-energy, 
with all integrations carried out by principal value, and 5^Beo is the corresponding contribution defined 
in Refs. [26, 27], due to squared imaginary parts. For 35-65 states, the above separation is not really 
essential, because S^Bqo is a numerically marginal contribution as compared to (see Ref. [26]), 
and thus bL{nS) ~ bL{nS) to a very good approximation. For P states under investigation here, the 
distinction (9), surprisingly, is already important (see Table 1). Final numerical values of the weighted 
difference of Beo for 5 states, and for individual P states, are summarized in Table 2. 



3. Status of the ^He+ Lamb Shift 

In the current section (we keep units with ?i = c = eo = 1), we would Uke to use the results 
described above, in order to infer the current theoretical status of the Lamb shift of 15 and 25 in 
the ^He+ ion. Before we start our actual discussion, however, we should remember that an ideal way 
to carry out a related calculation would involve a full-featured least-squares adjustment according to 
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Table 3. Contributions to ttie Lamb shifts of the lS'1/2 and 2S-1/2 states in ''He^. All equation numbers are 
connected to the contributions as listed in Ref. [7], unless indicated otherwise. SE= self-energy, VP= vacuum 
polarization, num. int.= numerical integration, rn = mass of orbiting particle, for the electron rn = nic, and 
M = nuclear mass. 



Order of contribution [rricc^} 


Equation in Ref. [7] 


£(15') [MHz] 


C{2S) [MHz] 


a(Zat\n\(Za)-'^] 

V / L\ / J 


Eq. (59) [part] 


146 724.762 


18 340.595 


a{Za)^ 


Eq. (59) [part] 


-40 796.296 


-4 725.621 




Eq. (59) fpart] 


16.295 


2.037 


a^{Zaf 


Eq. (59) [part] 


0.029 


0.004 


a{Za)'^ (muonic vac. pol.) 


Eq. (63) 


-0.081 


-0.010 


a{Za)'^ (hadronic vac. pol.) 


Eq. (65) 


-0.051 


-0.006 


a{Za)^ (SE+VP) 


Eq. (73) 


1 827.214 


228.402 


cx^ { Z(x) ^ (two one-loops) 


Eqs. (74-I-76-I-79-I-80) 


0.492 


0.061 


a^{Zaf (two-loop VP) 


Eq. (75) 


0.704 


0.088 


a^iZaf (two-loop SE) 


Eq. (81) 


-10.709(1) 


-1.339 


c^i^Zcef' (two-loop sum) 


sum of 3 above 


-9.513(1) 


-1.189 


a{Zaf\n^\{Za)-'^\ 


Eq. (84) [part] 


-198.172 


-24.771 


a{Zaf\Q[{Za)-^] 


Eq. (84) [part] 


123.906 


16.985 


a{Zaf'GsE(Za) 


Ref. [28] 


-82.542 


-10.621 


a{ZafGvp{Za) 


Ref. [29] 


-1.685 


-0.276 


a{ZafGwK{Za) 


Eq. (101) 


0.157 


0.020 


a^{Zafln^[{Za)-'^] 


Ref. [4] 


-1.153 


-0.144 


a^{Zafln^[{Za)-^] 


Ref. [6] 


-0.213 


0.037 


a2(Z«)«ln[(Za)-2] 


Ref. [8] 


2.666 


0.279 


a\Zaf 


Refs. [8, 14] 


-0.607(211) 


-0.064(26) 


{Zafm/M 


Eq. (136) 


17.786 


2.547 


{ZafmlM 


Eq. (144) 


-0.119 


-0.015 


{Zayin'^{Za)m/M 


Eq. (147) 


-0.010 


-0.001 


a{Zafm/M 


Eq. (151)+Eq. (46) of Ref. [30] 


-0.112 


-0.014 


a{Zaf\n\{Za)-'^]m/M 


Eq. (155) 


0.018 


0.002 


[Zafim/Mf 


Eq. (15) of Ref. [31] 


-0.053 


-0.007 


Z{Zaf{m/Mf 


approx., Eq. (152) 


0.019(19) 


0.002(2) 


Nucl. size [rel., 1.680(5) fm] 


Eq. (17) of this work 


70.865(422) 


8.860(53) 


Nucl. size [rel., 1.673(1) fm] 


Eq. (17) of this work 


70.275(84) 


8.786(11) 


Sum [1.680(5) fm] 




107693.112(472) 


13837.031(59) 


Sum [1.673(1)1111] 




i()7()l)2.r)22(228) 


l:!8:!().l)o7(29) 



Ref. [33], which includes all available data from relevant high-precision experiments (see [17]) and 
which, in principle, allows for a deduction of the nuclear charge radius. In order to infer an approximate 
theoretical prediction, though, one has to use a charge radius obtained from other sources, and we 
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Table 4. Contributions to tiie Lamb sliifts of the 2P-1/2 and IPzji states in '^He'''. As in Table 3, all equation 
numbers are connected to the contributions as listed in Ref. [7], imless indicated otherwise. The acronyms used 



for the corrections are also the same as in Table 3. 


Order of contribution [meC^] 


Equation in Ref. [7] 


-^(2^1/2) [MHz] 


'C(2P3/2) [MHz] 


a(ZQ)*ln[(Za)-^] 


Eq. (60) [part] 


0.000 


0.000 




Eq. (60) [part] 


-206.095 


200.725 




Eq. (60) [part] 


0.414 


-0.207 


oC'{Zaf 


Eq. (60) [part] 


-0.003 


0.002 


a(Zaf (muonic vac. pol.) 


Ref. [17] 


0.000 


0.000 


oi(Zoi)'^ (hadronic vac. pol.) 


Ref. [17] 


0.000 


0.000 


a{Z(xf (SE+VP) 


Eq. (73) 


0.000 


0.000 


Q^(Za) ' (two one-loops) 


Eqs. (74-h76-h79-h80) 


0.000 


0.000 


o?{Zaf (two-loop VP) 


Eq. (75) 


0.000 


0.000 


a^{Zaf (two-loop SE) 


Eq. (81) 


0.000 


0.000 


o?{ZoLf' (two-loop sum) 


sum of 3 above 


0.000 


0.000 


a{Zaf \n\{ZaY'^\ 


Eq. (86) [part] 


0.000 


0.000 


a{Zaf\x^\{ZaY'^\ 


Eq. (86) 


1.677 


0.944 


a{ZafG^Y.{Za) 


Ref. [28] 


-0.329 


-0.163 


a{ZafGwp{Za) 


Ref. [29] 


-0.022 


-0.005 


a{ZafGwi^{Za) 


Eq. (101) 


0.000 


0.000 


Q^{Zaf \vL\ZaY'^\ 


Ref. [4] 


0.000 


0.000 


O t r-w \Ri O r / r-w \ On 

a^{Zaf\n\Za) ^] 


Ref. [4] 


0.006 


0.006 


a^{Zaf\n[{Za)~'^] 


Ref. [8] 


0.001 


-0.001 


a^{Zaf 


Refs. [8, 16] 


-0.001 


-0.001 


{Zafm/M 


Eq. (136) 


-0.138 


-0.138 


{Zafm/M 


Eq. (145) 


0.007 


0.007 


{Zayin^{Za)m/M 


Eq. (147) 


0.000 


0.000 


ot{Zafm/M 


Eq. (151) 


0.000 


0.000 


a{Zaf\n^[{Za)-'^]m/M 


Eq. (155) 


0.000 


0.000 


{Za)^{rn/M)'^ 


E.g., Ref. [32] 


0.002 


-0.001 


Z{Zay'{rn/A'lf 


approx., Eq. (152) 


0.000 


0.000 


Nucl. size 


liq. (17) of this work 


O.OIK) 


0. ()()() 


Sum 




-2l)i.l8i 


2()i.i()8 



intend to follow this different route in the current work. 



We partly base our evaluation on Refs. [7,17,18] and choose a format as in Table 1 ofRef. [17]. In 
the evaluations described in Tables 3 and 4, the 2002 CODATA values of the fundamental constants [34] 
were used. 
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For the Lamb shift L, we use the imphcit definition [31, 32] 



E = m, [f{n,j) - 1] - [f{n,j) - 1]' + £ + E^, ■ 




(10) 



Here, E is the energy level of the bound two-body system under investigation, and /(n, j) is the dimen- 
sionless Dirac energy. E.g., we have /(1, 1/2) = /(15) = y^l - (Za)^, and /(2, 1/2) = f{2S) = 



y i (1 + y/l — {Za)^) for the 15 and 25' states, respectively. Furthermore, rrir is the reduced mass 
of the system, M is the nuclear mass, and Ehfs is the energy shift due to hyperfine effects, which are 
absent for the spinless ^Hc nucleus. 

In order to avoid confusion, we would like to include a few clarifying words regarding specific 
entries in Tables 3 and 4. In general, we have added the factor a to all scales for the contributions Usted 
in Tables 3 and 4. giving all contributions with an overall scaling of TOoC^, in contrast to a mcC?, which 
had been used in Ref. [17]. Regarding the contribution of order (Za)^ ln[(ZQ;)^^] in the first row of 
Table 3, it is worthwhile to note that this term represents the leading logarithm of the Lamb shift, given 



Here, rrii is the reduced mass of the system, given by nir = merriN / (we + ton), where is the mass 
of the nucleus. The reduced-mass dependence of the argument of the logarithm itself, ln[(Za)~^] — > 
ln[{Za)~^ nie/mr] = ln[(Za)~^] + ln[me/mr], is being included here and in Ref. [17] into the non- 
logarithmic term of order a{Za)'^. One might wonder why there is a theoretical uncertainty associated 
to this contribution at all. The reason is that the most accurate theoretical value for the quantity (1 1) is 
obtained by expressing it in terms of the 2002 CODATA Rydberg constant, which has a relative uncer- 
tainty of 6.6 X 10~^^, and the 2002 CODATA fine-structure constant, which has a relative uncertainity 
of 3.3 X 10~^^. The latter is responsible for the small theoretical uncertainty of the leading logarithmic 
contribution to the Lamb shift of the 15 level. The term of order aiZctf' contains both contributions 
from the self-energy and the vacuum polarization, as indicated by the explanatory note "SE+VR" 

The term of order aiZoLf ln^[(Za)~^] corresponds to the self-energy coefficient ^62, as given 
e.g. in Ref. [35], and the indicated term of order a{Z(xf' ln[(Za)^^] is the sum of a self-energy and a 
vacuum-polarization contribution in this order. Note that the latter distinction differs from the one used 
in Ref. [17], where the term of order a{Zoif W'\{Za)~'^\ denotes the sum of a double logarithmic, 
and a single logarithmic self-energy contribution, and the term of order a{Z(yf lii[(Za)~^] was re- 
served exclusively for the vacuum-polarization contribution in this order. The values of the self-energy 
remainder G^^{Za) for 5 and P states are listed in [28, 36]. The vacuum polarization remainder 
function Gypl^a) is taken from [29] and corresponds exclusively to the Uehling part of the one- 
loop vacuum polarization. It might be worthwhile to point out that at the current level of accuracy, it 
is entirely sufficient to consider the vacuum-polarization higher-order remainder for P states via the 
formula 




by 




(11) 



l^E^AnP,) = ^ f'-^Y {Ajo^inP,) + (Za) Aj,^ (nP,)} , 



(12) 



where the analytic coefficients read 



(13a) 



A^^{nPy2) 



4l7r n^-l 



A^^inPy^) 



777 



(13b) 



2304 n2 



768 n2 
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These have been obtained in Refs. [7, 16, 37] for general principal quantum number n. 

For the i362-term for S states of order a^{Za)^ In [(Za)^^], the result from Ref. [6] was used, 
which supersedes the estimate given in Eq. (101) of Ref. [7]. For the analytic Sgi-term of order 
a'^{Za) \n[{Za)~'^], the result given in Ref. [8] provides the most recent value. 

For the B^q coefficient corresponding to the nonlogarithmic term of order a^{Za)^ for the ground 
state, two mutually contradictory results of -61.6± 15 % (Ref. [25]) and -127±30 % (Ref. [14]) have 
been reported. The latter is based on an extrapolation of an all-order (in Za) numerical calculation. 
Note, however, that there is a known missing piece in the analytic result reported in Ref. [25], which 
is currently under study and which will need to be evaluated before final conclusions can be drawn. 
M. Bides [38] therefore suggested that a valid interim way of estimating the uncertainty of Bqq would 
consist in taking the arithmetic mean of these two results, and taking the half difference as an estimate 
for the theoretical uncertainty. This uncertainty would comprise all higher-order analytic terms, as it 
involves a comparison to a nonperturbative (in Za) calculation. For the 25' state, we can use the result 
for 15* and add the weighted difference listed in Table 2. For P states, the results reported in Sec. 2 
of this paper (see also [16]) provide enough information to eliminate all theoretical uncertainty at the 
current level of accuracy. 

Finally, let us remark that a term of order a'^{Zay hi^[{Za)^'^] could be estimated in principle 
on the basis of taking a "local" Lamb-shift potential that corresponds to the self-energy part of A50, 
namely, 



SV = 4a{Za)^ 



139 1 , 



(.4, 



taken as an input for a Dirac-^ induced correction to the one-loop self-energy. The result of this ap- 
proach could alternatively be used as an uncertainty estimate for all the higher-order terms. This pro- 
cedure leads to the estimate 



Br2{nS) = ±-7r 



139 1 , 

128 - 2 ^"(2) 



(15) 



For IS, this leads to a value of ±0.337 MHz for the higher-order two-loop remainder. However, the 
uncertainty due to the B72 -contribution is already contained in the uncertainty of the remainder term 
of order a^iZa)^, because in determining the uncertainty of Bqq, a comparison was made to a non- 
perturbative numerical calculation for higher nuclear charge numbers. The latter necessarily contains 
all contributions from the B72 term and aU higher-order remainders. The above result of ±0.337 MHz 
therefore likely overestimates the uncertainty due to the higher-order two-loop remainder and is men- 
tioned here only for illustrative purposes. 

Concerning radiative-recoil corrections, we note that the discrepancy between [39] and [40] con- 
cerning radiative insertions into the electron line in the a{Za)^m/M radiative recoil correction was 
resolved in [30] . The analytical result from that work was used. According to Eq. (46) of Ref. [30] , the 
coefficient multiplying the non- vacuum-polarization part of order correction of order a{Za)^m/M is 
-1. 32402796... /n^ 

The nuclear spin in ''He^ is different as compared to atomic hydrogen. The former is a spin- 
1/2-spin-zero system, whereas the latter is a spin- 1/2-spin- 1/2 system. Recoil corrections of first 
order in the mass ratio are unaffected by the different spin of the nucleus as compared to hydrogen. 
However, recoil terms of order Z{Za)^{m/M)'^ are nuclear spin-dependent. Without carrying out a 
detailed analysis, we approximately calculate the nuclear self-energy effects of order Z{Za)^ {m/M)^ 
by leaving out the Pauli form-factor correction from Eq. (153) of Ref. [7], which is certainly absent 
for a spinless nucleus, and we conservatively take the Dirac form factor contribution as an uncertainty, 
while we note that a more detailed analysis would be of interest and currently lacking in the literature. 



©2006 NRC Canada 



Jentschura and Haas 



9 



Finally, we add the nuclear-spin dependent correction listed in Eq. (15) of Ref . [31], 



AE = 



1 {Za)'^mc /"rric 



y 



(16) 



which is of second order in the mass ratio, for the spin-l/2-spin-zero system under investigation. This 
term is connected to the absence of the zitterbewegung term in the Breit Hamiltonian for a spinless 
nucleus. For P states, we also add terms of order (Za)'^ ( ' which do not depend on the zitterbe- 
wegung term and are given, e.g., in Ref. [32]. 

Concerning the nuclear-size correction, we would like to mention that a full integration of a nuclear 
potential with a fully relativistic wave function (e.g., within a hard-sphere approximation) turns out to 
be quite essential to obtain reliable values for this correction. We have carried out such an integration 
in the current investigation with the full Dirac wave functions and obtain results in agreement with 
Ref. [41]. Results obtained with two different values for the two different root-mean- square charge 
radii of 1.673(1) fm [42] and 1.680(5) fm [43] are given in Table 3. Note that the former charge radius 
of 1.673(1) fmhas been questioned (see e.g. [44,45]). The uncertainty due to the shape of the nuclear 
charge distribution can be estimated to be much smaller than the uncertainty due to the nuclear size, 
based on experience with highly charged ions [46]. For the ^He-nucleus, a spherically symmetric model 
is well justified (closed shell, spin / = 0). 

The nuclear-size correction AEfs{nS) and AEfs{nPj), for low nuclear charge numbers, can be 
approximated very well by the first few terms of an expansion in the two small parameters Za and 
m{r^y^^, with the result 



AEf,{nS) 



{ZaYm^R'' 



1 



{ZamR)+ 



+{Zaf 



In 



2 Z am R 



2j + *(n) 



227 
150 



2 

5 n 



10n2 J 



(17a) 



AEkinPi) =— ^ '— — 5j^i,2 . 

10 n-' ' 



(17b) 



Here, 7 = 0.577216 ... is Euler's constant, \l/(n) is the logarithmic derivative of the Gamma function, 
and R is the radius of the nucleus in a hard-sphere model, which is related to the root-mean-square 
radius (r^)^/^ by the following formula [see Eq. (7) of Ref. [41]], 



R 



2\l/2 



3/3 (r4) 1 



(18) 



For the 2Pi/2 state, we obtain an upward finite nuclear-size energy shift of 353 Hz which is barely 
significant on the kHz level (see Table 4). 

The results in Eqs. (17) and (18) have been obtained in the approximation of an infinitely heavy 
nucleus, and with exact Dirac wave functions for a point nucleus. Both of these approximations should 
be vaUd for ''He'*'. In addition, it should be noted that both the results given in Eqs. (17) and (18) 
are in excellent numerical agreement with a full numerical integration of the finite-size potential with 
Dirac wave functions. The linear correction term in R, i.e. the term — i (ZamR) in Eq. (17a), is 
a consequence of the exponential factor « exp^—Zamr/n) in the wave function, which should not 
be ignored in the evaluation of the finite-size effect, although this effect is primarily sensitive to the 
probability density at the origin (at the nucleus). Any further effects that influence the finite-size effect 
like nuclear polarization are here absorbed into the uncertainty of the nuclear radius (see, e.g., Ref. [47] 
for an illustrative discussion of some of the further aspects that are relevant to the finite-size effect). 
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4. Conclusions 

In Sec. 2, we summarize recent theoretical results for the higher-order two-loop binding corrections 
to the Lamb shift. These results are used, in Sec. 3, to infer updated values for the Lamb shift of low- 
lying states of the *He+ ion. Some of the analytic coefficients used in the evaluation are given in 
Eqs. (13) and (17). The analytic expansion of the nuclear finite-size correction (17) might be useful in 
other contexts. 

The recent progress in the field has allowed for an improvement of the theory beyond the limits 
set by the leading-order effects, and for some of the most accurate predictions in all of theoretical 
physics. In particular, we reemphasize that for P states, the results reported in Sec. 2 of this paper (see 
also [16]) provide sufficient information to eUminate all theoretical uncertainty at the kHz level for the 
Lamb shift in the ^He"*" ion. 
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